Search for charmonium and charmonium-like states in T(25) radiative decays 
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Using a sample of 158 million T(2S) events collected with the Belle detector, charmonium and 
charmonium-like states with even charge parity are searched for in T(25*) radiative decays. No 
significant \c.j or r\ c signal is observed and the following upper limits at 90% confidence level (C.L.) 
are obtained: B(T(2S) ->■ 7Xco) < 1.0 x 10" 4 , B(T(25) ->■ 7Xci) < 3.6 x 10" 8 , B(T(2S) -> ix^i) < 
1.5 x 10 -5 , and S(T(2S) — > 7?j c ) < 2.7 x 10~ 5 . No significant signal of any charmonium-like state is 
observed, and we obtain the limits B(T(25) -> 7^(3872)) x B(X(3872) ->• tt+tt" J/tji) < 0.8 x 10" 6 , 
B(T{2S) 7 X(3872)) x B(X(3872) -» tt+tt-tt J/i>) < 2.4 x 10" 6 , B(T(2S) -> 7 X(3915)) x 
B(X(3915) -> uJ/ip) < 2.8 x 10~ 6 , B(T(25) ->■ 7 F(4140)) x B(F(4140) -> < 1.2 x 10" 6 , 

and B(T(25) -> 7 X(4350)) x B(X(4350) < 1.3 x 10~ 6 at 90% C.L. 

PACS numbers: 14.40. Pq, 14.40.Rt, 13.20.Gd 



The data samples of the B factories have provided a 
wealth of experimental information on charmonium spec- 
troscopy [l|. Below open charm threshold agreement 
between experimental mass measurements and predic- 
tions based upon potential models was recently demon- 
strated with high accuracy for the h c [1, However, in 
the region above the open charm threshold, in addition 
to many conventional charmonium states, a number of 
charmonium-like states (the so-called "XYZ particles") 
have been discovered with unusual properties. These may 
include exotic states, such as quark-gluon hybrids, meson 
molecules, and multi-quark states (lj. Many of these new 
states are established in a single production mechanism 
or in a single decay mode only. To better understand 
them, it is necessary to search for such states in more 
production processes and/or decay modes. States with 
can be studied via initial state radiation 
(ISR) with the large T(4S') data samples at BaBar or 
Belle, or via e + e~ collisions directly at the peak energy 
at, for example, BESIII. For charge-parity-even charmo- 
nium states, radiative decays of the narrow T states be- 
low the open bottom threshold can be examined. 

The production rates of the P-wave spin-triplet 
XcJ (J=0, 1, 2) and S- wave spin-singlet r) c states in 
T(15) radiative decays have been calculated by Gao et 
ai; the rates in T(25') decays are estimated to be at the 
same level Q. However, there are no such calculations 
or estimations for "XYZ particles" due to the limited 
knowledge of their nature. 



In this paper, with the world largest data sample taken 
at the T(2S) peak, we report a search for the Xcj, ?7c, 
X(3872) [|, X(3915) @, and F(4140) in T(2S) ra- 
diative decays, extending our previous work on the T(15) 
sample 8]. In addition, the new structure X(4350) 0], 
which was observed as a 3.2 standard deviation (a) signal 
in 77 — ¥ <f>J /ip is also searched for. As any charmonium 
state above il>(2S) is expected to have a larger branch- 
ing fraction for the El/Ml transition to 'il^(2S) than to 
J/ip [l(|, we also search for states decaying into jtp(2S). 

The data used in this analysis include a 24.7 fb _1 
data sample collected at the T(25) peak and a 1.7 fb _1 
data sample collected at y/s = 9.993 GeV (off-resonance 
data) with the Belle detector o pera ting at the KEKB 
asymmetric-energy e + e~ collider [12j . The number of 
the T(2S*) events is determined by counting the hadronic 
events in the data taken at the T(25*) peak after sub- 
tracting the scaled continuum background from the data 
sample collected at </s — 9.993 GeV. The selection 
criteria for hadronic events are validated with the off- 
resonance data by comparing the measured R value (R = 
*(e+e-->hadronsU ^ C LEO's result fil . The number 

er(e + e —S-fi^fj, ) 1 !— 

of T(25) events is determined to be (158 ±4) x 10 6 , with 
the error dominated by the MC simulation of the T(2S) 
decay dynamics using pythia (lif . 

Well measured charged tracks and photon candidates 
are first selected. For a charged track, the impact pa- 
rameters perpendicular to and along the beam direction 
with respect to the interaction point (IP) are required 
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to be less than 0.5 cm and 4 cm, respectively, and the 
transverse momentum should exceed 0.1 GeV/c in the 
laboratory frame. Information from different detector 
subsystems is combined to form a likelihood Li for each 
particle species [15j . A track with TZk — Ck+c > 0-6 
is identified as a kaon, while a track with IZk < 0.4 is 
treated as a pion. With this selection, the kaon (pion) 
identification efficiency is about 90% (96%), while 5% 
(6%) of kaons (pions) are misidentified as pions (kaons). 
For electron identification, the likelihood ratio is defined 
as TZ e = c , where C e and C x are the likelihoods for 
electron and non-electron, respectively, determined using 
the ratio of the energy deposited in the electromagnetic 
calorimeter (ECL) to the momentum measured in the sili- 
con vertex detector and central drift chamber (CDC), the 
shower shape in the ECL, the matching between the po- 
sition of charged track trajectory and the cluster position 
in the ECL, the hit information from the aerogel thresh- 
old Cherenkov counters and the dE / dx measurements in 
the CDC [lH . For muon identification, the likelihood ra- 
tio is defined as = c +C f ' +CK , where £ M , £ w and Ck 
are the likelihoods for muon, pion and kaon hypotheses, 
respectively, based on the matching quality and pene- 
tration depth of associated hits in the iron flux return 
(KLM) [ljj . A good neutral cluster is reconstructed as a 
photon if its ECL shower does not match the extrapola- 
tion of any charged track and its energy is greater than 
40 MeV. In the e + e~ center-of-mass (CM.) frame, the 
photon candidate with the maximum energy is taken to 
be the T(25') radiative decay photon (denoted as tr), 
and its energy is required to be greater than 3.5 GeV. A 
3.5 GeV photon energy corresponds to a particle of mass 
5.5 GeV/c 2 produced in T(2S*) radiative decays. 

We reconstruct J /ijj signals from e + e~ or candi- 
dates. In order to reduce the effect of bremsstrahlung 
or final-state radiation, photons detected in the ECL 
within 0.05 radians of the original e + or e~ direction 
are included in the calculation of the e + /e~ momen- 
tum. For the lepton pair used to reconstruct J/ip, 
at least one track should have 7Z e > 0.95 while the 
other should satisfy TZ e > 0.05 in the e + e~ mode; or 
one track should have > 0.95 (in the XcJ analy- 
sis, the other track should have associated hits in the 
KLM detector that agree with the extrapolated trajec- 
tory of a charged track provided by the drift chamber) 
in the mode. The lepton pair identification ef- 

ficiency is about 97% for J/tp — > e + e~ and 87% for 
J r /ip — > fi + fi~ . In order to improve the J/ip momen- 
tum resolution, a mass-constrainted fit is then performed 
for J/ip signals in the jJ/ip, ir + ir~ J/ip, ir + ir~ ir° J / ip , 
and (pJ/ip modes. Different modes have similar J /ip 
mass resolutions. The J/ip signal region is defined as 
\Mi+i- — m J/Tp\ < 30 MeV/c 2 (« 2.5er), where mj/^ is 
the nominal mass of J/ip- The J/ip mass sidebands are 
defined as 2.959 GeV/c 2 < M e+e - < 3.019 GeV/c 2 and 
3.175 GeV/c 2 < M l+l - < 3.235 GeV/c 2 , and are twice 
as wide as the signal region. For the jip(2S) channel, the 
-0(25*) is reconstructed from the ir + ir~J/ip final state, 



with a mass constrained to the nominal ip(2S) mass to 
improve its momentum resolution. To estimate the dif- 
ference in the ip(2S) mass resolution between MC simula- 
tion and data, the process e + e~ — > -fisRip{2S) is selected 
as a reference sample, and the mass resolution is found 
to be 3.0 ±0.1 MeV/c 2 from data, and 2.6 MeV/c 2 from 
MC simulation. The difference in the mass resolution is 
included when extracting the signal yields in the analyses 
below. 

We search for the \ c j in the "fJ/ip mode. The en- 
ergy deposited by the \cj photon (denoted as 7/, since 
its energy is much lower than that of 7r) is required to 
be greater than 150 MeV to reduce the large background 
from mis-reconstructed photons, and the total number of 
photons is required to be exactly two to suppress multi- 
photon backgrounds. The angle between the 7r and 7/ 
should be larger than 18° to remove the background from 
split-off fake photons. To remove the ISR background 
e + e~ — > 7isrV'(2S') — > 7isr7XcJ, where a photon is miss- 
ing, we require the square of the "mass recoiling against 
the 7/ and J/ip" (M 2 CC = (P e+e - - P f ) 2 , here P e+e - is 
the 4-momentum of the e + e~ collision system, and Pf 
is the sum of the 4-momenta of the observed final state 
particles) to be within -0.5 GeV 2 /c 4 and 0.5 GeV 2 /c 4 . 
This M 2 ec requirement is effective since this background 
has at least two missing photons and M 2 0C (7/ J/ip) tends 
to be large. Bhabha and dimuon background events with 
final-state radiative photons are further suppressed by re- 
moving events in which a photon is detected within a 18° 
cone around each charged track direction. 

The mode shows a clear J/tp signal, while the 

e + e~ mode has some residual radiative Bhabha back- 
ground. Figure [T] shows the jiJ/ip invariant mass distri- 
bution together with the background estimated from the 
J I ip mass sidebands (normalized to the width of the J / ip 
signal range) for the combined e + e~ and modes 
after the above selection criteria are applied. Some ISR 
backgrounds with a correctly reconstructed J/tp remain 
in the data. No XcJ signal is observed. 

A simultaneous fit to the signal region is performed 
with Breit-Wigner (BW) functions convolved with Gaus- 
sian resolution functions for the resonances and a second- 
order polynomial background term. The width of the 
Gaussian resolution function is fixed at 7.9 MeV/c 2 , 
which is obtained by increasing the MC-simulated value 
by 10% to account for the difference between data and 
MC simulation. The masses and widths of the \c.J res- 
onances are fixed to their PDG values [l9j]. In the si- 
multaneous fit, the ratio of the yields in the two J/tp 
decay channels is fixed to where is the J/ip de- 
cay branching fraction for the e + e~ mode or fi + /i~ mode 
reported by the PDG [19(, and Ei is the MC-determined 
efficiency for this mode. The upper limit on the number 
(n up ) of signal events at the 90% C.L. is calculated by 

/ C{x)dx 

solving the equation ° +OQ = 0.9, where x is the 

J C(x)dx 

number of signal events, and C{x) is the likelihood func- 
tion depending on x from the fit to the data. The values 
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of n up are found to be 2.8, 3.1 and 7.6 for the Xco, Xd an d 
Xc2i respectively, when requiring the signal yields to be 
non-negative in the fit. We do not observe any structure 
at high masses, where excited Xcj states are expected. 
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FIG. 1: The "fiJ/tp invariant mass distribution in the T(2S) 
data sample. There is no Xco, Xa> ° r Xc2 signal observed. The 
solid curve is the best fit, the dashed curve is the background, 
and the shaded histogram is from the normalized J/ip mass 
sidebands. The signal yield is required to be non-negative in 
the fit. 

To search for a possible excited charmonium state 
in the jiip(2S) final state, a J/ip candidate and two 
oppositely charged pion candidates are reconstructed. 
The ip(2S) signal region is defined as 3.67 GeV/c 2 < 
M^+„-j/,i, < 3.70 GeV/c 2 , and the ip(2S) mass side- 

< 
< 



bands are defined as 3.63 GeV/c 2 < M„+ n -j/^ 
3.66 GeV/c 2 and 3.71 GeV/c 2 < M n+w -j/^ 
3.74 GeV/c 2 . To suppress backgrounds with miscon- 
structed photons, we require the energy of the 7; to be 
higher than 75 MeV. To suppress the ISR background 
e + e~ — > 7isrV'(2S') — > 7isr7t + 7t ~J/ip> we require the 
square of the mass recoiling against the ~yi and ip(2S) to 
be within -0.5 GeV 2 /c 4 and 1.5 GeV 2 /c 4 since M 2 oc for 
the ISR background tends to be shifted towards negative 
values. 

The r )ii]j{2S) invariant mass distribution after the 
above selection is shown in Fig. [5] There is no signifi- 
cant signal. However, a few events accumulate around 
3.82 GeV/c 2 .where the 7^(25) decays of the Xco(2P) 
andr7 C 2(l-D) [lCj are expected. A fit between 3.75 GeV/c 2 
and 3.90 GeV/c 2 with a Gaussian to parameterize the 
signal shape yields a mass of (3.824 ± 0.002) GeV/c 2 and 
a signal yield of 5.5 ± 2.7 events corresponding to a sta- 
tistical significance of 1.8o\ The signal significance is 
determined by comparing the value of — 2 ln(Lo/^max) 
from the fit, with values from fits to 10,000 pseudo- 
experiments. Here Lq and L max are the likelihoods of 
the fits without and with the signal, respectively. The up- 
per limit on the product branching fraction B(T(2S) — > 
7A) x B(X 7^(25)) < 1.3 x 10~ 5 at the 90% C.L. is 
determined following the procedure described below. 

To search for the r] c signal in T(2S) radiative decays, 
we reconstruct rj c candidates from the K s K + ir~ + c.c, 
ir+n-K+K-, 2(K+K-), 2(7r+7r"), and 3(7r+7r-) modes. 
Well measured charged tracks should be identified as pi- 
ons or kaons, and the number of charged tracks is six 
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FIG. 2: The 7z - i/>(2S l ) invariant mass distribution. The open 
histogram is from the ip{2S) signal mass region, the shaded 
histogram is from the normalized ip(2S) mass sidebands. In 
the inset, the solid curve is the best fit between 3.75 GeV/c 2 
and 3.90 GeV/c 2 , and the dashed curve is a fit with only a 
second-order polynomial to describe the background. 



for the 3(7r + 7r~) final state and four for the other final 
states. In the KgK + tt~ + c.c. mode, Kg candidates are 
reconstructed from tt + it~ pairs with an invariant mass 
M„+„- within 30 MeV/c 2 of the K s nominal mass. A 
Kg candidate should have a displaced vertex and flight 
direction consistent with a K s originating from the IP; 
the same selection method is used in Ref. [llj]. Events 
with leptons misidcntified as pions in the ir + ir~ K + K~ 
and 2(7r + 7r~) modes are removed by requiring lZ e < 0.9 
and TZft < 0.9 for the pion candidates. The value of Af 2 c 
for the hadronic daughters of the rj c candidate is required 
to be within -1 GeV 2 /c 4 and 1 GeV 2 /c 4 . 

After the selection described above, Fig. [3] shows the 
combined mass distribution of the hadronic final states 
for the five r\ c decay modes. The large J/ tp signal is due 
to the ISR process e + e~ — > 7isr<//'0> while the accumu- 
lation of events within the r\ c mass region is small. The 
shaded histogram in Fig. [3] is the same distribution for 
the off-resonance data and is not normalized. 
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FIG. 3: The mass distribution for a sum of the five r\ c de- 
cay modes. The solid curve is a sum of the corresponding 
functions obtained from a simultaneous fit to all the r) c de- 
cay modes, and the dashed curve is a sum of the background 
functions from the fit. The shaded histogram is a sum of the 
off-resonance events (not normalized). The J/ip signal is pro- 
duced via ISR rather than from a radiative decay of an T(nS) 
resonance. 
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A simultaneous fit is performed to the five final states. 
The ratios of the r\ c (J/ip) yields in all the channels are 
fixed to BiSi, where each B% is the r\ c (J /if)) decay branch- 
ing fraction for the i-th mode reported by the PDG [l9j . 
and Ei is the MC-dctermincd efficiency for this mode. 
The fit function contains a BW function convolved with 
a Gaussian resolution function (its resolution is fixed to 
7.9 MeV/c 2 from MC simulation) describing the rj c sig- 
nal shape, another Gaussian function describing the J/ip 
signal shape, and a second-order polynomial describing 
the background shape. The mass and width of the BW 
function are fixed to the PDG values for the r\ c . The 
results of the fit are shown in Fig. [3j where the solid 
curve is the sum of all the fit functions, and the dashed 
curve is the sum of the background functions. The fit 
yields 1 4 ± 20 r\ c signal events corresponding to an upper 
limit n up of 44 at the 90% C.L. In addition, we obtain 
370±15 J/tp signal events from the fit (in agreement with 
338 ± 16 expected from jisnJ/ip production according to 
MC simulation), giving a mass of 3098.1 ± 0.7 MeV/c 2 , 
which is consistent with the PDG value (l9j . 

The selection criteria for T(2S) -> 7rX(3872), 
A (3872) —> ir + it~ J/tp are similar to those used for ISR 
ir + n~J/tp events in T(4S f ) data [20]. We require that 
one J/tp candidate be reconstructed, two well-identified 
7r's have an invariant mass greater than 0.35 GeV/c 2 , 
and that M? cc (ir + ir~ J/ip) be within the range between 
-1 GeV 2 /c 4 and 1 GeV 2 /c 4 . To suppress the ISR 
7r + 7r~ J/tp background, we require that the polar angle of 
the 7r candidate satisfy |cos6*| < 0.9 in the e + e~ CM. 
frame. Except for a few residual ISR produced tp(2S) 
signal events, only a small number of events appear in 
the TT + TT~J/ip invariant mass distribution, as shown in 
Fig. 0Ja). There is no accumulation of events in the 
A (3872) mass region. Fitting using a signal shape from 
the MC sample and a first-order polynomial function as 
the background shape, the upper limit n up for the num- 
ber of signal events is determined to be 3.6 at the 90% 
C.L. 

We also search for the A(3872) and A(3915) in the 
ir + ir~n J/tp mode. We select 7r + , 7r~, and J/ip can- 
didates in the A(3872) —> n + ir~J/tp mode (with the 
requirement on the ir + ir~ invariant mass greater than 
0.35 GeV/c 2 removed) and a n° candidate from a pair 
of photons with invariant mass within 10 MeV/c 2 of the 
7T° nominal mass. Here the ir° mass resolution is about 
4 MeV/c 2 from MC simulation. Figure SJb) shows the 
7T + 7r~7r° J/tp invariant mass distribution, where the open 
histogram is the MC expectation for the A (3872) sig- 
nal plotted with an arbitrary normalization. Using the 
same fit method as in A(3872) —> it + it~ J/tp, we deter- 
mine n up for the number of A (3872) signal events to be 
4.2 at the 90% C.L. Figure QJc) shows the scatter plot 
of m(jr + TT~TT°J/tp) versus m(7r + 7r~7r°) from data, where 
the region indicated by the ellipse corresponds to the 
±3<7 mass regions of m(7r + 7r _ 7r° J/tp) and m(7r + 7r~7r°) 
from the A(3915) — !• ujj/tp decay. There is one event 
with m(n + TT-TT°J/tp) at 3.923 GeV/c 2 and m(7r+7r-7r°) 
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FIG. 4: (a) Distribution of the tt+tt" J/tp 
for T(25) -> 7r7t + 7t~ J/tp candidates, (b) 
the 7r + 7r _ 7r° J/tp invariant mass for T(2S) — > 7R7r + 7r _ 7r° J/tp 
candidates. (c) Scatter plots of m(7r + 7r _ 7r° J/tp) versus 
m(7T + 7r~7r°), where the region indicated by the ellipse cor- 
responds to the ±3ct mass regions of m(ir + iv~ iv J/tp) and 
m(-K + -K~-K°) from the X(3915) —¥ uij/tp decay. Points with 
error bars are data, open histograms are the MC expectation 
for the X(3872) signal (arbitrary normalization). The peak 
at 3.686 GeV/c 2 in (a) is due to tp(2S) production via ISR. 



at 0.790 GeV/c 2 from T(2S) data, as shown in the el- 
lipse. Assuming that the number of background events 
is zero, the upper limit n up for the number of A(3915) 
signal events is 4.4 at the 90% C.L. 

We search for the F(4140) and the A(4350) in the 
4>J/tp mode. The selection criteria are very similar to 
those in the analysis of AT(3872) — > ir + ir~ J/tp described 
above and the <p is reconstructed from a K + K~ pair. 
According to MC simulation, the (p signal region is de- 
fined as 1.01 GeV/c 2 < M K+K - < 1.03 GeV/c 2 . The 
number of well measured charged tracks is required to be 
exactly four. After applying all of the above event selec- 
tion criteria, there is no clear J/tp or <p signal. Nor are 
there candidate events in the y(4140) or A(4350) mass 
regions. The upper limits on the number of V(4140) and 
A(4350) signal events are both 2.3 at the 90% C.L. 
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Several sources of systematic uncertainties are con- 
sidered. The uncertainty due to particle identification 
efficiency is 2.4%-3.4% and depends on the hnal state 
particles. The uncertainty in the tracking efficiency for 
tracks with angles and momenta characteristic of signal 
events is about 0.35% per track, and is additive. The 
photon reconstruction contributes an additional 2.0% per 
photon. Errors on the branching fractions of the inter- 
mediate states are taken from the PDG [f9]; they are 
6.9% for the Xco mode, 4.5% for the Xd mode, 4.2% 
for the Xc2 mode, 1.7% for the jip(2S) mode, 17% for 
the rj c mode, 1.0% for the A(3872) mode, 1.3% for the 
A(3915) mode, and 1.6% for the </>J/V> mode. By using a 
phase space distribution and including possible interme- 
diate resonant states, the largest difference of efficiency 
is determined to be 2.1% for the r\ c decay modes. The 
difference in the overall efficiency for a flat angular distri- 
bution of radiative photons and a l±cos 2 9 distribution is 
less than 3.0%. Therefore, we quote an additional error 
of 5.0% due to the limited knowledge of the decay dy- 
namics for all the states studied, except for the Xco mode 
and rj c mode, which are known to follow a 1 + cos 2 
distribution. According to MC simulation, the trigger 
efficiency is 89% for the Xcj mode, rather high for other 
modes (> 99%); we take a 3.0% error for the XcJ mode 
and 1.0% error for other modes as a conservative esti- 
mate of the corresponding uncertainties. With the pure 
e + e~ — » 7isr?/'(2S'), V'(2S') — > Tr + ir~J/ip or J/ij}rj{-^- 77) 
samples obtained from Belle data, the uncertainty due to 
the recoil mass squared requirement is 1.0% for the chan- 
nels with a single photon and 4.7% for channels with two 
photons. By changing the order of the background poly- 
nomial, the range of the fit, and the values of the masses 
and widths of the resonances, uncertainties on the XcJ 
and r\ c signal yields are estimated to be 1.1% and 16%, 
respectively. In the T(2S) — > JrXcJ mode, the uncer- 
tainty associated with the requirement on the number of 
photons is 2.0% after applying a correction factor of 0.94 
to the MC efficiency, which is determined from a study 
of a very pure T(2S) — > event sample. In the 

rjc K s K + n~ + c.c. mode, the uncertainty in the K s 
selection efficiency is determined by a study on a large 
sample of high momentum Kg — > iy + tt~ decays; the ef- 
ficiency difference between data and MC simulation is 
less than 4.9% [2l|. Finally, the uncertainty on the total 
number of T(2S) events is 2.3%. Assuming that all of 
these systematic error sources are independent, we add 
them in quadrature to obtain a total systematic error as 
shown in Table HI 



Since there is no evidence for signals in the modes stud- 
ied, we determine upper limits on the branching fractions 
of T(25 l ) radiative decays. Table fl] lists the upper limits 
n up for the number of signal events, detection efficien- 
cies, systematic errors, and final results for the upper 
limits on the branching fractions. In order to calculate 
conservative upper limits on these branching fractions, 
the efficiencies are lowered by a factor of 1 — cr sys in the 
calculation. 

TABLE I: Summary of the limits on T(25) radiative decays to 
charmonium and charmonium-like states R. Here n up is the 
upper limit on the number of signal events, e is the efficiency 
with the secondary decay branching fractions excluded and 
trigger efficiency included, cr sys is the total systematic error, 
and B(T(2S) -> 7_R) up (Br) is the upper limit at the 90% 
C.L. on the decay branching fraction in the charmonium state 
case, and on the product branching fraction in the case of a 
charmonium-like state. 



State (R) 


n up 


e{%) 


0- sys (%) 


Br 




XcO 


2.8 


14.2 


10.9 


1.0 x 10" 


-4 


Xcl 


3.1 


14.8 


10.8 


3.6 x 10" 


■6 


Xc2 


7.6 


15.2 


10.7 


1.5 x 10" 


-5 


r]c 


44 


26.3 


24 


2.7 x 10~ 


-5 


X(3872) n+ir-J/tf) 


3.6 


27.3 


7.4 


0.8 x 10- 


■6 


X(3872) -> ir + ir-ir°J/iP 4.2 


10.3 


9.6 


2.4 x 10- 


-6 


X(3915) -> ujj/i} 


4.4 


10.5 


9.6 


2.8 x 10- 


■6 


F(4140) -> (j)J/ii 


2.3 


22.3 


7.4 


1.2 x 10- 


■6 


X(4350) ->■ (pJ/ip 


2.3 


21.0 


7.4 


1.3 x 10" 


■6 



To summarize, we find no significant signals for the XcJ 
or rjc, as well as for the A(3872), A(3915), F(4140), or 
A(4350) in T(25') radiative decays. The results obtained 
on the Xcj and i] c production rates are consistent with 
the theoretical predictions of [|J . 
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